Deterministic generation of Gaussian pure state in quasi-local dissipative system 
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This paper shows that an arbitrary Gaussian pure state can be deterministically generated in 
a dissipative open system that has quasi-local interactions between the subsystems and couples 
to surrounding environment in a local manner. A quasi-local interaction, which means that the 
interaction occurs among only a few subsystems, is a crucial requirement for practical engineering 
of a dissipative system. The key idea is that, first, an auxiliary system having local interaction 
with environment is prepared, and then that auxiliary system is coupled to the underlying target 
system via a set of two-body Hamiltonians in such a way that a desired pure state is generated. 
Moreover, we show that, even with a simple single- mode auxiliary system, deterministic generation 
of an arbitrary approximate Gaussian cluster state is possible, by devising an appropriate switching 
scheme. Lastly, we discuss in a specific example how much a dissipation-induced pure Gaussian 
state can be perturbed by decoherence and parameter error. 

PACS numbers: 03.65.Yz, 42.50.Dv 



I. INTRODUCTION 

Preparing a desired pure state is a crucial task in quan- 
tum information technologies. However, in a realistic sit- 
uation any quantum system unavoidably interacts with 
surrouding environment and generally obeys the follow- 
ing Markovian master equation: 



d 
dt 



p(t)=-i[H,p(t)] 



m f 1 1 \ 

where H is the system Hamiltonian and Li is the coupling 
operator representing the interaction between the system 
and the i-th environment channel. In usual, the state p(t) 
in Eq. (pQ) dissipatively evolves in time towards a mixed 
state and never recover its purity. However, recently it 
has been shown in several papers [H-Q that, by engineer- 
ing suitable pairs of the system operators H and L^ the 
state p(t) can be uniquely moved to a pure steady state. 
This means that a desired pure state can be deterministi- 
cally generated without specific initialization of the sys- 
tem. We actually find some applications of this environ- 
ment engineering approach to entangled state generation 
[8-16] , and further, some advanced quantum information 
processing such as quantum computation [17j, memories 
[HI, and distillation 

Here we mention about the quasi-local interaction. An 
operator is called quasi-local if it acts only on a given 
site of the system and its neighborhood; particularly in 
this paper we call the interaction quasi-local if it occurs 
among only two nodes. In the case of a quantum oscil- 
lator network whose i-th node is a single-mode oscilla- 
tor of variable xi = (qi,Pi), an example of a quasi-local 
coupling is L = + <72- However, to generate a desired 
pure state such as a highly entangled state via the above- 
mentioned dissipation-based method, the system is of- 
ten required to have non quasi-local (global) interaction 
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FIG. 1: (a) Example of a quantum oscillator network having 
non quasi- local interaction with environment. The number i 
denotes the i-th oscillator, (b) An extended system that yields 
the same pure Gaussian steady state as that produced by the 
system (a). The auxiliary system couples to the environment 
locally, while it couples to the target system in a quasi-local 
manner through a certain Hamiltinian. 



with environment. Let us consider a 7- nodes quantum 
oscillator network depicted in Fig. 1 (a); this network 
has internal couplings among the nodes xs-xq and X4-X5 
through a Hamiltonian H, and further, it has two global 
couplings of the nodes (#1, #2, #3, £4) an d (#5, #6? #7) to 
the environments, the operators of which are for instance 
L\ = qi + q2 + q3 + q± and L 2 = P5 + P6 + P7- From a 
practical viewpoint, clearly, this kind of global interac- 
tion is hard to implement. Actually there have been some 
proposals to engineer a dissipative system having quasi- 
local interaction with environment, which yet determin- 
istically produces a useful pure state. For instance, it was 
shown in [To| that a qubit 2D cluster state can be gen- 
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erated via a quasi- local dissipative process, and also Ref. 
[liSj demonstrated the case of a uniformly distributed en- 
tangled state. Moreover, Ticozzi and Viola showed in 
[20| that, in the general finite dimensional case, it is al- 
ways possible to construct a dissipative system generat- 
ing a desired pure state as well as satisfying a certain 
quasi- locality constraint. However, in the infinite dimen- 
sional case, only a few specific results have been reported 

The result obtained in [6] is that, for any pure Gaus- 
sian state, we can always engineer an open Gaussian sys- 
tem whose unique steady state is identical to that pure 
state. However, as mentioned above, this open Gaussian 
system often has to have global interaction with envi- 
ronment, such as the system shown in Fig. 1 (a). In 
this paper, nonetheless, we show that an arbitrary pure 
Gaussian state can be dissipatively generated in a cer- 
tain extended system having only local interaction with 
environment. The idea is described as follows: First, we 
couple the underlying target system, via a set of two- 
body Hamiltonians, to a certain auxiliary system having 
local interactions with the environment. That is, the ex- 
tended open system composed of the target and the aux- 
iliary systems contains only quasi- local interactions; see 
Fig. 1 (b). Then, based on the result of [6], we will show 
that a Gaussian pure state generated in the original sys- 
tem of interest, which can include global interactions with 
environment, is also generated in the above-mentioned 
extended system. This configuration is, actually, a gen- 
eralization of the result of [8[, where dissipation-induced 
generation of a two-mode squeezed state in atomic en- 
sembles was demonstrated. 

In addition, we present a simple schematic that only 
utilizes a single-mode auxiliary system. More specifically, 
even in the case where the above-mentioned result re- 
quires us to have a multi-mode auxiliary system, it will 
be shown that only a single-mode auxiliary system can 
serve to achieve the same goal. The key idea is the use 
of a switching scheme of the auxiliary system; that is, 
in the case of the system depicted in Fig. 1, a single- 
mode auxiliary system first couples to the subsystems 
xi, . . . , £4, and then, it switches to couple to the subsys- 
tems X5,X6,X7. A particularly important result is that, 
with this switching scheme, an arbitrary approximate 
Gaussian cluster state |2ll-[25| can be generated dissipa- 
tively. This is a generalization of the switching scheme 
proposed by Li et ai in [12], which demonstrated dis- 
sipative generation of several types of four-mode cluster 
states in an actual physical setup. 

Lastly, we examine how much a specific dissipation- 
induced pure state is robust against some perturbation; 
the system is a pair of two atomic ensembles that dissi- 
patively generates a two- mode squeezed state [8[, where 
here damping decoherence and parameter uncertainty are 
further taken into account. Then we find the optimal 
system parameters that maximize the entanglement be- 
tween the atomic ensembles. A point is that these opti- 
mal parameters are out of the range where the auxiliary 



cavity mode can be adiabatically eliminated; that is, this 
is an example where the extended system is really robust 
against some perturbation compared to only the target 
system obtained through the adiabatic elimination. 

Notation. Let X = (Xij) be a matrix whose entry Xij 
is an operator on a Hilbert space or a complex number. 
We then define the following matrices: 

X T = (X Jl ), xt = 

Sft(X) and $s(X) denote the real and imaginary parts 
of X. For a vector of operators x, we define the anti- 
commutator of x by 

{x,x T } = xx T + (xx T ) T . 

II. GAUSSIAN DISSIPATIVE SYSTEM 

In this section, we describe a general Gaussian dissi- 
pative system in terms of a quantum stochastic differen- 
tial equation (QSDE) and present the condition for the 
steady state of this system to be pure @. 

The time evolution of a general open quantum system 
is generated by the unitary operator U(t) subjected to 
the following QSDE (in Ito form) [26|: 

m 

du ^) = { ( L i d Mty - LtdAiit)) 
i=i 

Y^-L*L t +iH\dt}u(t), (2) 

with ?7(0) = I (we have not included the scattering 
term). Ai(t) is the annihilation process on the i-th envi- 
ronment vacuum field; this satisfies the quantum Ito rule, 
e.g., dAi(t)dAj(t)* = Sijdt. Let p be the initial state of 
the system and |0) be the vacuum state of the whole en- 
vironment field. The master equation (pQ) is obtained by 
differentiating p(t) = U{t)(p <8> |0)(0|)l7(t)* and tracing 
out the field state. 

Next, let (qi,Pi) be the canonical conjugate pair of 
the i-th quantum oscillator, which satisfies the canonical 
commutation relation (CCR) qiPj —pjqi = iSij. Defining 
the vector of observables x = [q±, . . . ,q n ,pi, . . . ,p n ] T ? we 
can write the CCR as 

xx {xx ) zXl n , S n j q , 

where I n denotes the n x n identity matrix (we often 
drop the subscript n). A general linear open system can 
be characterized by the following Hamiltonian H and the 
coupling operators L = [Li, . . . , L m ] T : 

H = -x T Gx, L = Cx, 
2 

where G = G T £ R 2 ™ x2 ™ and C e C mx2n > From Eq ^ 
the QSDE ofx(t) = [U "(*)* qi U (£) , U '(*)* piU (t) ,• • • ] T is 
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given by 



A. The extended dissipative system 



dx(t) = Ax(t)dt + BdW(t), 



(3) 



where A = £ n (G + C T S m C/2) and B = E n C T with 
C = V / 2[3?(iC) T ,^(zC) T ] T . Also, we have defined the 
field operators Q< = (A { + A*)/y/2, Pi = (A { - A?)/V2i, 
and W= [Qi,...,Q m ,Pi,...,P m ] T . 

A Gaussian system can be fully characterized by only 
the mean vector (x) = [(gi), (pi), • ■ • ] T and the co- 
variance matrix V = ({Ax, Ax T }} /2 with Ax = x — 
(x). Note that a Gaussian state is pure if and only if 
1/ Y // 2 2n det(V) = 1. The importance of the linear system 
(j3j) is in the fact that, if the initial state of this system is 
Gaussian, then it preserves the Gaussianity of the state 
for all time. In particular, the covariance matrix V(t) 
obeys the following Lyapunov differential equation: 

lv(t) = AV(t) + V(t)A T + ±BB T . (4) 

Hence, the covariance matrix of the steady state is given 
by the solution of the algebraic Lyapunov equation 



AV + VA 1 



1 



BB l 



0. 



(5) 



This equation has a unique solution iff A is Hurwitz; i.e., 
all eigenvalues of A are in the open left half complex 
plane. 

Now we are interested in a linear Gaussian system 
whose steady state is uniquely pure. It has been shown in 
[H that an arbitrary pure Gaussian state can be dissipa- 
tively generated, if the system matrices G and C can be 
freely chosen. In particular, the following result is useful 
in this paper: 

Theorem 1 [6]: Suppose that Eq. (j5j) has a unique 
solution V . Then, this is the covariance matrix of a pure 
Gaussian state if and only if the following matrix equa- 
tions are satisfied: 



V + -E n J C T = 0, Yi n GV + V(E n G) T = 0. (6) 



III. PURE GAUSSIAN STATE GENERATION 
VIA QUASI-LOCAL DISSIPATIVE 
ENVIRONMENT 

To generate a certain desirable pure Gaussian state 
dissipatively, the system is often required to have global 
interaction with the environment; that is, for instance, 
the coupling operator must be of the form L = q\ + 
Q2 + Q3 + • • which corresponds to C = (1, 1,1,.. .). In 
this section, we show that such a global interaction be- 
tween the system and the environment can be always 
avoided by constructing an auxiliary system that quasi- 
locally couples to the target system and locally couples 
to the environment. 



In addition to the underlying target system con- 
sisting of n quantum oscillators, we consider an aux- 
iliary system consisting of m oscillators. Let x = 
[gi, . . .,q n ,Pi, • • • ,Pn] T and x = [<?i, . . .,q m ,pi, • • • ,Pm] T 
be the vectors of canonical conjugate pairs of each sys- 
tem. We assume that the decoherence of the target sys- 
tem is negligible, i.e., the system operators are 

H = -x T Gx, L = 0, 
2 

where G = G T G ]R 2nx2n . The auxiliary system has no 
self-Hamiltonian and dissipates through typical damping 
channels, hence the system operators are given by 

H = 0. L = v^a, 

where a = [ai, . . . , a m ] T is a vector of annihilation oper- 
ators, i.e., di = (qi + ipi)/y/2. The two systems couple 
via the following interaction Hamiltonian: 



i(d)Cx — x T C^a), 



(7) 



with C G C mx2n . Note that H{ nt describes a quasi- local 
interaction between the two systems, because it can al- 
ways be decomposed into a sum of two-body Hamiltoni- 
ans as follows: 



where x^ = [qi,Pi] T and x^ = [^,Pi] T . We will see 
later that this quasi-local Hamiltonian ([7]) has the same 
effect as the coupling operator of the system of originally 
interest, i.e., L = Cx, when choosing G and C so that 
they satisfy the condition of Theorem 1. Defining C = 
v / 2[5R(^C) T ,^(iC) T ] T , we can rewrite Eq. © as 



-(x T C T x + x T Cx). 



[x(t) T ,x(t) T ] T 
d£(t) = A£(t)dt + BdW(t), (8) 



Then, the overall system vector £(t) 
obeys the following linear QSDE: 



where 



^rnC —tvl2m/^ 



B 







(9) 



The covariance matrix V(t) = ({Af(t), A£(t) T }}/2 
evolves in time through Eq. (jU), and its steady solution 
is, if it exists, obtained by solving Eq. (J5j). To avoid 
confusion, we again present the same equation 



AV + VA T + \BB T 



0, 



(10) 



where A and B are now given in Eq. (j9]). 

Recall here that we are originally interested in the sys- 
tem with Hamiltonian H = x T Gx/2 and the coupling 
operator L = Cx, which thus obeys the linear dynamics 

dx(t) = A lX {t)dt + BxdWQ), (11) 

where A x = £ n (G + C T E m C/2) and B 1 = £ n C T . Note 
that the dissipation channel L = Cx can be global. Espe- 
cially, its steady covariance matrix is our concern, which 
is subjected to 

A 1 V 1 + V 1 Aj + ±B 1 Bj = 0. (12) 

The following theorem states that, if the system ([TT]) has 
a pure steady state, then the extended system © can 
dissipatively produce the same pure Gaussian state. 

Theorem 2: Suppose Eq. ([T2]) has a unique solution 
V\ that corresponds to a pure Gaussian state. Then, 
Eq. ([TP]) has a unique solution V = diag(Vi, l2™/2). 

Proof. First, we prove that Eq. ([TO]) has a unique solu- 
tion, which is equivalent to that A is Hurwitz. For this 
purpose, let us set A T rj = At?; i.e., rj G C 2 ^ n+m ^ is an 
eigenvector of A T and A G C is the corresponding eigen- 
value. Then, multiplying Eq. ([TP]) by from the left and 
by 7? from the right, we have 5ft (A) = —rfBB T r]/ (4rfVrj). 
This further becomes Sft(A) = —^772 772 / (477"^ V77), where we 
have defined 77 = [rjj, 77 J] T with 771 G C 2n and 772 G C 2m . 
Now let us assume 7/2=0; then, from A T r] = A77 we 
obtain 

GE^=Aryi, CE^i = 0. (13) 

Moreover, noting the assumption that Eq. (fl2j) has a 
unique solution 

we find that 771 satisfies 

1 r°° 

v\Vim = 2 1 \\C^e A ^ m \\ 2 dt. 

But this takes zero due to Eq. ([T3]) and the Cayley- 
Hamilton theorem. This conclusion contradicts to V± > 
0, hence we have 772 7^ 0; this further leads to Sft(A) < 0, 
implying that A is Hurwitz. 

Next let us prove that V = diag(Vi, hm/ty is the solu- 
tion of Eq. ([TO]) . When representing V in a block matrix 
form V = [Vi, V12; V12, V2], these entries satisfy 

S„GT/i + Fi(S„G) T + E„C T ^ + yi 2 (S m (7) T = 0, 
^(S m C) T + Z n C T V 2 + E n GV 12 - kV 12 /2 = 0, 
S m CVi2 + (X m C) T - nV 2 + K 7 2m /2 = 0. (14) 



Now the assumption is that V\ corresponds to a unique 
pure steady state, hence from Theorem 1 we have 

^(S m (7) T + ^S„(7 T = 0, S n G^ + y(S n G) T = 0. 

(15) 

Note that the former condition comes from the fact that 
Theorem 1 yields V^iCY - Z n $(iC) T /2 = and 
V^iiCy + £ n J?(7C) T /2. Then, we find that the set 
of matrices V\2 = 0, V2 = ^2m/2, and V\ satisfying 
Eq. is the solution to Eq. (fH|) . Because Eq. ([H]) 
has a unique solution as shown in the former part of this 
proof, V = diag(Vi,l2m/2) is the unique solution. □ 

This theorem states that, once we find a suitable pair 
of matrices G and C such that the system ([TT]) dissi- 
patively and uniquely generates a desired pure Gaussian 
state, the same goal can be achieved by alternatively con- 
structing the extended system dSJ. Note again that the 
system ([TT]) can couple to the environment globally, while 
the extended system flHJ) locally couples to the environ- 
ment through the operator Li = ^fnai and its internal 
modes quasi-locally couple with each other through the 
Hamiltonian Hi nt . That is, the target system stabiliz- 
ing a desired pure state can be realized as a subsystem 
of an extended system having local coupling to the en- 
vironment and quasi-local internal couplings among the 
nodes. 

Remark 1: As an auxiliary system, we usually take 
a system with very fast modes that can be adiabatically 
eliminated. In our case, even when the assumption of 
Theorem 2 does not hold, by taking k sufficiently large, 
the auxiliary system rapidly converges to the vacuum and 
the mode x(t) can be adiabatically eliminated. Then, 
the target system is approximated b y th e system whose 
coupling operator is L = ICxj^Jn [27j. In this sense, 
Theorem 2 implies that we can treat the system as if the 
auxiliary modes were heavily damped as long as the pure 
steady state condition is satisfied. 



B. Generation of two-mode squeezed state 

We here consider a two-mode Gaussian system stud- 
ied in [8]. Although for this system the requirement of 
quasi-locality is already satisfied, this example clearly il- 
lustrates our idea. 

The physical setup is depicted in Fig. [2j the target sys- 
tem is two atomic ensembles trapped in a ring-type cav- 
ity, while the auxiliary system corresponds to this optical 
cavity composed of two propagating modes. These two 
systems interact with each other via external pulse lasers 
with Rabi frequencies fl u . and Q Si , (i = 1,2). Here we 
assume that the number of atoms in each ensemble is suf- 
ficiently large; then, the collective spin component of the 
atomic ensemble can be approximated by an annihilation 
operator, and consequently the interaction Hamiltonian 
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FIG. 2: Two atomic ensembles trapped in a two- modes op- 
tical cavity. Q. Ui and Q Si are the Rabi frequencies of the z-th 
external laser fields. 



Now, for the matrix C given by Eq. (p~6|) and G = 
0, the Lyapunov equation (fl2j) has the following unique 
solution: 



cosh(20 - sinh(2£) 
- sinh(2£) cosh(2£) 







cosh(20 sinh(20 
sinh(2f) cosh(2f) 



where £ = tanh _1 (r). This is the covariance matrix of 
a pure two- mode squeezed state. Therefore, from The- 
orem 2, the Lyapunov equation (fTD|) has a unique solu- 
tion V = diag(Vi, I4/2); that is, the pair of atomic en- 
sembles acquires the two-mode squeezed state at steady 
state, while the auxiliary cavity mode becomes a trivial 
coherent state. Equivalently, the whole four- mode atom- 
cavity system generates the same atomic steady state as 
that generated in the two atomic ensembles system with 
Hamiltonian H = and the coupling operators 



is given by 

H int = {aliy/flhPuiCLi + \fN~ 2 f3 S2 a* 2 ) + H.c.} 

+ {a* 2 (^/m/3 Sl al + y/N2/3 U2 a 2 ) + H.c.} , 

where a$ and a$ are the annihilation operators of the i-th 
atomic ensemble and the i-th cavity mode, respectively. 
Ni denotes the number of atoms of the i-th ensemble, 
and also 



a, 



2A,, 



9 8* 



2A, 



1,2, 



where g 9 and A. denote the coupling strength and the 
detuning, respectively. If we set the parameters as Ni = 
N,g Ui = 9si = g (i = 1, 2), and A n = A s = A, then the 
interaction Hamiltonian can be simply written as 



^± [ K(^ Ml ai + ^ S2 a*) + H.c.} 
+ {a* 2 (Q Sl al + n u2 a 2 ) + H.c.} ] . 



n S2 =n>o 



Further, let us set Rabi frequencies as Q Ul 
and Q U2 = Q 3l = r^, where r G [0, 1) is a parameter. 
Then, the interaction Hamiltonian is of the form ([7j) with 



V2 L 



1 r 
r 1 



i 

-ir 



-ir 



(16) 



where ji = \^NgQ/2A. For 87 Rb atoms, the sponta- 
neous emission of each atom is negligible and thus we 
can set L = 0, for typical values of the parameters, /i, hi w 
100 kHz and r = 0.8 [![. Also, the number of atoms is 
enough large so that the self Hamiltonian of the atomic 
ensembles is assumed to be 0, i.e., H = or equivalently 
G = 0. 

The coupling operator of the cavity is given by L = 
yfhza, where k is the damping rate, and we assume that 
the detuning of the cavity is 0, i.e., H = 0. 



Li = /i(ai + ra 2 ), L 2 = fi(a 2 + ra\). 

Another physical realization of this purely dissipative sys- 
tem was proposed in [HI . 



IV. DISSIPATION-INDUCED PURE CLUSTER 
STATE WITH SINGLE ENVIRONMENT 
CHANNEL 

As shown in the previous section, an arbitrary pure 
Gaussian state can be generated in the target system, by 
introducing a certain m-mode auxiliary system that lo- 
cally couples to m environment channels. From engineer- 
ing viewpoint, it is clearly convenient if we can achieve 
the goal with the auxiliary system having a small num- 
ber of modes m. In particular, let us consider the simple 
case H = 0; i.e., the system Hamiltonian H is negligible 
compared to the interaction Hamiltonian, as in the case 
of Sec. IIII B[ In this case, we need the condition m = n, 
which is understood from Theorem 1 together with the 
fact that, for a covariance matrix V corresponding to 
a pure Gaussian state, V + iS n /2 has an n-dimensional 
kernel [28|. However, this requirement is very demanding 
especially when n is large. 

In this section, we present a general solution to the 
above-posed problem. That is, we prove that, even for a 
general n-mode system with H = 0, only a single-mode 
auxiliary system coupling to a single environment chan- 
nel introduces a dissipative mechanism that drives the 
system state to an arbitrary approximate Gaussian clus- 
ter state (2ll-[25|. This result is significant in the sense 
that a most simple dissipative system deterministically 
generates a most useful quantum state from the quan- 
tum information viewpoint. The key idea is the use of 
a switching scheme of the interaction Hamiltonian be- 
tween the target and the auxiliary systems, which is a 
generalization of the idea proposed by Li et ai in [l2| . 
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A. The switching scheme 

As introduced above, in this section we consider a most 
simple system; that is, the system's Hamiltonian is neg- 
ligible (H = 0) and the auxiliary system is single-mode 
(m = 1), in addition to the assumptions L = H = and 
L = y/Ra. We represent the interaction Hamiltonian (O 
in the following form: 

n 

Him =^Y1 { h *Mj a 3 + Pjtf) + HX '}' ( 17 ) 

i=i 

where aj (j = l,...,n) denotes the j-th annihilation 
operator of the target system. otj,(3j G C, and /iGl are 
parameters. 

First let us consider the case n = 1 and set the pa- 
rameters as a± = 1 and /?i = r G [0,1). Then, the 
corresponding C matrix in Eq. ([7]) is given by C = 
/i[l+r, z(l — r)]/\/2 and Eq. ([T2]) has the following unique 
solution: 



-2€ o 







=2? 



(18) 



where £ = tanh _1 (r). This is the covariance matrix of a 
pure squeezed state, hence from Theorem 2 we find that 
Eq. ([TO]) has a unique solution V = diag(Vi,l2/2). As 
a result, when n = 1, the target system acquires a pure 
squeezed state while the auxiliary state becomes vacuum, 
through the interaction Hamiltonian ([17]) . 

Next we consider the case of n > 1. In this case, 
as mentioned above, Eq. ([T2]) does not have a unique 
solution, hence Theorem 2 cannot be directly applied. 
Now let us consider a unitary matrix U GC nxn and the 
unitary coordinate transformation a' = Ua of the vector 
of system variables a = [a±, . . . ,a n ] T . In terms of the 
quadratures q[ = (a - + a-*)/\/2 and p\ = (a - - a'*)/>/2i, 
we find that a/ = [q f \p f ] is the symplectic and orthogonal 
transformation of x as follows: 



5a;, 5 



'St(U) -$*(Uj 



Note that, since 5 is symplectic, x' satisfies the CCR: 
xx (x x ) — zS^jfiS — — 

We here explain the idea of switching. At the k-th 
switching stage, the parameters (a^ k \f3j k ^) of the inter- 
action Hamiltonian ([TT]) are chosen as follows: 



a 



(k) 



u kj , p. 



(k) 



rU* kj1 re [0,1), 



(20) 



where Ukj is the (k,j) element of the unitary matrix in- 
troduced above. With this choice, the interaction Hamil- 
tonian ([T7]) is written as 



H, 



(k) 



int -/iK(4 + r4*) + H.c.}. (21) 



This is no more than the Hamiltonian discussed in 
the case n = 1. Hence, in the k-th. switching stage 
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system 
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Vacuum 



(qi,Pi}\ — Q— q — o — o — { 



-o- 

2nd 



3rd 



Target system 



-O- 

4th 



Squeezed 



FIG. 3: The switching scheme. 



the k-th mode (q' k1 p' k ) deterministically changes to a 
pure squeezed state with covariance matrix ([15]) . Note 
that during the k-th stage, the other system variables 
(QeiP'e)) ^ 7^ ^ no ^ change at all. Therefore, ap- 
plying the interaction Hamiltonian ([2T]) repeatedly by 
changing k = 1,2, . . . , n, the corresponding system vari- 
able (q' k1 p' k ) gets squeezed in this order, as schematically 
shown in Fig. [3] In particular if the initial state is the 
ground state with respect to x' (and thus x as well), 
all the off-diagonal elements of the covariance matrix re- 
main and then the system's covariance matrix V{ = 
({Ax' J Ax' T }}/2 becomes V{ = diag(e- 2 ^/ n , e 2 ^I n )/2. 
As a result, the steady state of the target system of the 
mode x = S T x f is a pure Gaussian state with covariance 
matrix 











(22) 



This is a symplectic transformed pure squeezed state; we 
will see later that this state can represent any approx- 
imate Gaussian cluster state by appropriately choosing 
the unitary matrix U, or equivalently the switching pa- 



(19) rameters (af\pf ] 



)• 



B. Generation of CV cluster state 

A Gaussian cluster state is an entangled state of great 
importance, particularly in one-way quantum computa- 
tion [22|, [25|| . Hence, this state should be a target that is 
dissipatively generated with the scheme presented in this 
paper, particularly with the switching scheme. Thus we 
here show how to chose the unitary matrix U so that the 
covariance matrix ([22]) represents a given target Gaussian 
cluster state. 

Here is the definition of a Gaussian cluster state: Let 



A = A G 



be the adjacency matrix representing 



the graph structure of a cluster state of interest; i.e., 
the element of A represents the weight of the cou- 
pling between the i-th and the j-th nodes of the network. 
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Then, the approximate Gaussian cluster state is defined 
as a state satisfying 



lim cov(p — Aq) = 0, 



(23) 



where a G R is a certain parameter contained in the state; 
it often corresponds to a squeezing parameter. Also we 
have defined cov(x) = ({Ax, Ax T })/2. It is of course 
impossible in reality to take the limit a —> oo, hence 
we call the state with finite a the approximate Gaussian 
cluster state. 

Now we describe a relation that connects a given adja- 
cency matrix A and the unitary matrix U characterizing 
the covariance matrix ([22]) . As will be mentioned later 
in Remark 2 that one such relation has already been ob- 
tained in [23|, [lH particularly for the aim of constructing 
a concrete optical process corresponding to U ', but we 
here provide an alternative method in rather an abstract 
way: 



Proposition 1. Let A = A T G 



be the adjacency 



matrix of a given Gaussian cluster state and define N — 
— {il n + A). Then the polar decomposition N = RU 
yields a unitary matrix U and a real matrix R. The 
Gaussian state with covariance matrix ([22]) characterized 
by this unitary U then satisfies 



cov(p-Aq) = -(I n +A 2 )e~ 2 ^ 



(24) 



hence it is an approximate Gaussian cluster state con- 
verging to the idel one in the limit of £ — >> 00. 



Proof. First, R is real, because R = ViVF = I n + A 2 G 
R nxn . Next, for the Gaussian state with covariance ma- 
trix (l22l) we have 



cov(p — Aq) 



e" 2 «/„ 
e 2 «/„ 

1 

2 





-A 


s 





where Fi = $l(U)A+$t(U) and F 2 = %(U)A-$t(U). Not- 
ing that RU = N and R is real, we have 3?(C7) = 
and 3?([7) = -iT 1 . Hence, we have F ± = -R~ l (I n + A 2 ) 
and F2 = 0. Furthermore, noting that I n = UW = 
R- 1 NN*R- T , we have = AWt = 7 n + A 2 , thus 

F^F 1 = (I n + A 2 )i?- T i?- 1 (/ n + A 2 ) =I n + A 2 . As a 
result, we have cov(p — Aq) = (J n + A 2 )e~ 2 ^ /2. □ 

The merit of this result is that the unitary matrix [/ is 
straightforwardly constructed from a given A compared 
to the result of [23|, [25| , although in this case U does not 
have a clear correspondence to some optical realizations. 
However, Eq. ([20]) clarifies how to physically implement 
the interaction Hamiltonian ([T7]) . 

Remark 2: The approximate Gaussian cluster state 
with covariance matrix ([22]) can be deterministically gen- 
erated on optical fields, by the following method |23l. l25j. 
We first prepare n independent and identical squeezed 



light fields, and mix them via some passive optical de- 
vices such as a beam splitter and a phase shifter in a 
specific order determined from the adjacency matrix A. 
The collection of these transformations is totally repre- 
sented by a unitary matrix; if we denote that unitary 
matrix as U\ the covariance matrix of the output fields 
is identical to Eq. ([22]) . The biggest difference of this 
setup to the presented scheme is that here we consider 
a state generated in a matter while the above-mentioned 
optical state lives in a flying light field; the former can 
be later manipulated or stored, while the latter is suited 
for propagating quantum information. A specific relation 
between these two regimes was discussed in [29| . 

Remark 3: In Proposition 1 the polar decomposition 
is used to find the appropriate unitary matrix U, but it 
is clear from the proof that only a certain decomposition 
of the form Af = RU with U unitary and R real gives the 
same relation ([24]) . Here, we see that, besides the polar 
decomposition, the Gram Schmidt procedure also serves 
as a convenient method to obtain such a decomposition. 
Let U be a unitary matrix whose row vectors are obtained 
from the Gram Schmidt procedure of the row vectors of 
N. Then, it is immediate to have the relation A" = RU 
with a real lower triangular matrix R, which is called the 
RQ decomposition. Hence, the Gram Schmidt procedure 
also yields the unitary matrix satisfying the condition in 
Proposition 1. 



C. Physical realization 

We can implement the presented switching scheme for 
n atomic ensembles (target system) trapped in a single- 
mode optical cavity (auxiliary system), in a similar con- 
figuration studied in Sec. IIIIB1 The physical setup in the 
case n = 4 is depicted in Fig. HJ Let aj be the annihila- 
tion operator approximating the collective spin compo- 
nent of the j-th atomic ensemble and a\ be the annihila- 
tion operator of the cavity mode. Then, the interaction 
Hamiltonian is given by [12| 



Ng 



2A 



[aUSlu^aj + n Sj e^a*) + H.c. 



where <p 9 G [0,27r) is the laser phase. The switching 
scheme shown in Sec. IIVI suggests us, at the k-th switch- 
ing stage, to chose the parameters as 

n uj e**"t =SlU kj , n Sj e 1 ^ =rnU£ p re [0,1), (25) 

where Q > is a parameter. The unitary matrix (Ukj) 
is determined from the target Gaussian cluster state. As 
proven in Sec. IIV A[ the whole state of the atomic en- 
sembles deterministically reaches this target, if they are 
all initially set to the ground states. 
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1 cavity field 
tt Ul ,Q Sl (auxiliary) 

(94, P4) 



QhPi 



4 atomic 
ensembles 
(target) 



J ^2 



(92,02) 




Pa) 



FIG. 4: Atomic ensembles trapped in a single- mode opti- 
cal cavity. Q Ui and Q Si are the Rabi frequencies of the z-th 
external laser fields. 



D. Example 

Let us consider the problem of dissipatively generating 
a four-nodes square cluster state depicted in Fig. [5j in 
the atomic ensembles system discussed in the previous 
subsection. The connecting edges between the nodes are 
equally weighted, and thus the adjacency matrix A of 
this graph state is given by 



A 



1 r 
0011 
1100 
1100 



We follow Proposition 1 to determine a unitary matrix U 
characterizing the switching law ([20]) . In particular, here 
we utilize the Gram Schmidt procedure; see Remark 3. 
To make the calculation simple, we orthogonalize the fol- 
lowing matrix: 



N' 



"1 -1 





- 




~i 


—i 





" 





1 


-1 


TV = 








i 


—i 


1 








i 





1 


1 


_o 


1 


_ 




1 


1 


i 






Then, U is given by 



U = 



-»[i,-i,o,o]/V?" 
-»[o,o,i,-i]/V5 

2,2]/ VlO 
[2,2,M]M0 . 



(26) 



Actually with this choice we have 



cov(p — Aq) 



■3 2 0- 

2 3 

3 2 

2 3. 



(27) 



FIG. 5: 

state. 



Graph structure of the four-nodes square cluster 



thus the state well approximates the target square cluster 
state when large £ is taken. 

The switching law ([25]) with the unitary matrix ([26]) 
clarifies how to choose the laser parameters as shown 
below. At the first switching stage, they are determined 
from the first row vector of U as 



Qui 



1 

= 72 
0, i 



ft, z = l,2, 
3,4, 



7T, 4>s 



1 



-7T. 



Through the interaction Hamiltonian with these parame- 
ters, the CCR pair {q'i,p'i) gets squeezed in the long time 
limit. Second, we switch the parameters and set: 



Sis- 1 

Q — _£i — __o 



3,4, 



2*' 



bs z = ^7T, 



1 

= 2 n > 



3 

= 2 W - 



Then the second node (q^P^) gets squeezed. By re- 
peating a similar procedure, each CCR pair 
i = 1,...,4 becomes a squeezed state with covariance 
matrix ([T8j) . As a result, if each atomic ensemble is in the 
ground state at the initial time, the whole state changes 
to a pure Gaussian state with covariance matrix ([22]) . 
implying that the system reaches the target approximate 
Gaussian cluster state satisfying Eq. ([27]) . 

All the above results, including the switching law of the 
parameters of the interaction Hamiltonian, have been ob- 
tained in ; here we have shown the same result from 
a general standpoint. That is, from the general theory 
developed in Sees. II V Al and IIVB[ now we know that any 
approximate Gaussian cluster state can be deterministi- 
cally generated in atomic ensembles trapped in a single- 
mode cavity. Moreover, as demonstrated here, the appro- 
priate switching law can be systematically constructed, 
once the target cluster state is specified. 



V. PERTURBATION TO PURE STEADY STATE 

In this section, we reconsider the two atomic ensem- 
bles discussed in Sec. lIIIBl with taking into account some 
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specific perturbations added to the system. Actually, an 
atomic ensemble constructed in a cavity often loses co- 
herence due to spontaneous emission. In addition to this 
kind of quantum effect, in practice any system contains 
some parameter uncertainties, which can bring serious 
loss of coherence as well when aiming to dissipatively 
generate a pure state. In the case of an atomic ensem- 
ble, the number of atoms is usually never determined 
exactly. We take these two typical perturbations and 
evaluate how much the steady state is affected by these 
losses. In particular, we find the optimal squeezing level 
and the cavity damping rate that maximize the entan- 
glement. Although these investigations do not straight- 
forwardly provide a new insight on the quasi-locality dis- 
cussed throughout this paper, the result will clarify a 
merit of enlarging the system rather than focusing only 
on the target system obtained by adiabatic elimination. 



A. Decoherence effect 

First we assume that the atomic ensembles are sub- 
jected to loss of atomic coherence, in which case the 
corresponding coupling operator is represented by L = 
\/7[ a i? a 2] T with 7 the decoherence rate. The coefficient 
matrix of the driving term of the system, given in Eq. ((9]), 
is then changed to 



A' = 



£ 2 C T " 




'-lh/2 0" 


_E 2 C -Kh/2 








= A + AA, 



where C is defined via Eq. (fTO|) . The eigenvalues of A are 
given by 



A. 



16/i 2 (l -r 2 ). 



A notable point is that, when large squeezing is intro- 
duced, i.e., r w 1, the eigenvalue A + approaches to zero, 
even in the case 7 <C /i, n. This loss of stability of the 
system implies that a desirable convergence of the state is 
prevented. Thus, the perturbation AA is not negligible 
particularly when aiming to generate a large entangled 
state. Now the coefficient matrix of the diffusion term of 
the system is B' = diag(y / 7/4, ^JTlIa) and the Lyapunov 
equation A'V + V A ,T + B'B' T /2 = has the following 
explicit solution: 



Vi = diag 



rc 



-1/2 

-c rc 



1/2 



rc+ 1/2 c 
c rc + 1/2 



where c = 4rft/(ft + 7) {^7 + 4(1 — r 2 )}. Here, for sim- 
plicity we have replaced ft//i and 7//1 by k and 7, re- 
spectively. The entanglement of this two mode Gaussian 
state can be quantified by the logarithmic negativity [30| 
and is given by 

En = max{0, — log(2z/)}, 

_ ^7 2 + 4«(1 - r) 2 + 7(V + 4(1 - r 2 )) 
V ~ 2(ft + 7)^7 + 4(1 -r 2 )) ' 




FIG. 6: Logarithmic negativity En versus the squeezing level 
£ = tanh -1 (r) and the cavity damping rate for the cases 
(a) 7 = 0, (b) 7 = 0.01, and (c) (7, e) = (0.01, \flX). 



Fig.[6]shows En versus the squeezing level £ = tanh _1 (r) 
and the cavity damping rate ft, for the cases where the de- 
coherence rate 7 takes (a) 7 = and (b) 7 = 0.01. In the 
former case (a), since there is no decohering effect, the 
steady state is pure with maximal entanglement En = 2£ 
without respect to K] i.e., large squeezing directly means 
large entanglement. However, once the system is sub- 
jected to the decoherence, the steady state entanglement 
drastically degrades, even for very small decoherence rate 
7 = 0.01, as depicted in Fig.[6](b). In particular, it is ap- 
parent from the figure that large squeezing brings about 
large loss of entanglement, which is yet consistent to the 
fact that a large squeezed state is usually very fragile to 
decoherence. It is also reasonable that now n affects En 
in two ways; large n induces rapid leaking of the pho- 
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tons through the cavity light field, while small k means 
that the cavity mode interacts with the atomic ensembles 
many times, implying many emissions of the photons into 
the cavity field. Consequently, there exists an optimal set 
of the parameters n and r = tanh(£) that maximizes En' 

where d = [-J 6 + 5 7 4 - 2 7 2 + ( 7 2 + 1) >/-7 4 (7 2 ~ 4 )] 1/3 - 
Particularly in the case 7 = 0.01, these optimal values 
are given by r* w 0.97 and w 0.52, and then takes 
about 2.91. 

Remark 4: It should be pointed out that the above 
optimal parameters are out of the range where the aux- 
iliary cavity mode can be adiabatically eliminated; that 
is, as discussed in Remark 1 in Sec. IIII A\ the auxiliary 
system can be adiabatically eliminated only when k is 
sufficiently large, but clearly in this case En goes down 
to zero. In this sense, the system studied here provides 
an example where the extended system is really robust 
against decoherence compared to only the target system 
obtained through adiabatic elimination. 

B. Parameter uncertainty 

Since the condition imposing the system to have a pure 
steady state is described by a set of algebraic equations, 
it is easily violated by some parameter changes. In 
the case of the atomic system under consideration, the 
numbers of atoms of each ensembles must be exactly the 
same, i.e., N\ = A^, but it is fairly unrealistic. Hence let 
us here examine how much the difference of N\ and N2 
affects the enta nglemen t of the steady state. We partic- 
ularly set e = yj N2 /Ni = \/LT; in usual, the number of 
trapped atoms is of order 10 6 , and such relatively large 
uncertainty (10%) can actually happen. Fig.[6](c) shows 
the logarithmic negativity En, where the decoherence 
due to the spontaneous emission discussed in the pre- 
vious subsection is additionally taken into account. As 
easily expected, further degradation of the entanglement 
is observed, and in almost all ranges of the parameters £ 



and k the state is not anymore entangled. Nevertheless, 
surprisingly, it is not a uniform degradation; actually 
as in the previous case, by engineering the system with 
optimal parameters we obtain a steady state 

that still has relatively large entanglement of En = 2.41, 
which is only 17% loss of entanglement compared to the 
ideal value En = 2.91 where there is no uncertainty 
(i.e., e = 1). In other words, by constructing the system 
with these parameters (fi£,rj), we can guarantee the 
entangled state with at least En = 2.41, against the 
uncertainty of the difference of the number of atoms 
up to 10%. This robustness property indicates possible 
effectiveness of the dissipation-based method for state 
preparation even in a realistic situation, as long as the 
system parameters are appropriately determined. 



VI. CONCLUSION 

The main results of this paper are twofold: First, we 
have shown that an arbitrary Gaussian pure state can be 
deterministically generated via local dissipative environ- 
ment, by constructing an appropriate auxiliary system. 
Second, we have shown that, even when only a single- 
mode auxiliary system is available, a well-tuned switch- 
ing scheme allows us to stabilize any approximate Gaus- 
sian cluster state in a dissipative way. The former is a 
generalization of the scheme proposed in [8| that yields a 
dissipation-induced two- mode squeezed state, while the 
latter is that of [HI where deterministic generation of 
several four-mode cluster states was demonstrated. 

The essential mechanism for bringing quasi-locality to 
the system considered in this paper is that, for a Gaus- 
sian system, any interaction Hamiltonian is always a sum 
of two-body (hence quasi-local) Hamiltonians. This im- 
plies that, as long as this kind of interaction Hamiltonian 
is taken when constructing an auxiliary system, any non- 
Gaussian (target) system couples to environment quasi- 
locally. This would be an interesting approach to ex- 
plore a general method of constructing a desired quasi- 
local dissipative environment for general non- Gaussian 
systems. 
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